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ization requires that the fields have nonminimal couplings to the spacetime Ricci 
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inflation. The nonminimal couplings induce a nontrivial field-space manifold in the 
Einstein frame, and they also yield an effective potential in the Einstein frame with 
nontrivial curvature. The ridges or bumps in the Einstein-frame potential can lead 
to primodial non-Gaussianities of observable magnitude. We develop a covariant 
formalism to study perturbations in such models and calculate the primordial bis- 
pectrum. As in previous studies of non-Gaussianities in multifield models, our results 
for the bispectrum depend sensitively on the fields' initial conditions. 
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I. INTRODUCTION 

Inflationary cosmology remains the leading account of the very early universe, consistent 
with high-precision measurements of the cosmic microwave background radiation (CMB) 
[H 12]. A longstanding challenge, however, has been to realize successful early- universe 
inflation within a well-motivated model from high-energy particle physics. 

Realistic models of high-energy physics routinely include multiple scalar fields [3l H]. 
Unlike single-field models, multifield models generically produce entropy (or isocurvature) 
perturbations. The entropy perturbations, in turn, can cause the gauge-invariant curvature 
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perturbation, (, to evolve even on the longest length-scales, after modes have been stretched 
beyond the Hubble radius during inflation [3HT2] . Understanding the coupling and evolution 
of entropy perturbations in multifield models is therefore critical for studying features in the 
predicted power spectrum, such as non-Gaussianities, that are absent in simple single-field 
models. (For reviews see [UHlG].) 

Recent reviews of primordial non-Gaussianities have emphasized four criteria, at least 
one of which must be satisfied as a necessary (but not sufficient) condition for observable 
power spectra to deviate from predictions of single-field models. These criteria include 
[16]: (1) multiple fields; (2) noncanonical kinetic terms; (3) violation of slow-roll; or (4) 
an initial quantum state for fluctuations different than the usual Bunch-Davies vacuum. As 
we demonstrate here, the first three of these criteria are generically satisfied by models that 
include multiple scalar fields with nonminimal couplings to the spacetime Ricci curvature 
scalar. 

Nonminimal couplings arise in the action as necessary renormalization counterterms for 
scalar fields in curved spacetime [T7H22] . In many models the nonminimal coupling strength, 
^, grows without bound under renormalization-group flow [20]. In such models, if the non- 
minimal couplings are ^ ~ C>{1) at low energies, they will rise to ^ ^ 1 at the energy scales 
of early-universe inflation. We therefore expect realistic models of inflation to incorporate 
multiple scalar fields, each with a large nonminimal coupling. 

Upon performing a conformal transformation to the Einstein frame — in which the grav- 
itational portion of the action assumes canonical Einstein-Hilbert form — the nonminimal 
couplings induce a field-space manifold that is not conformal to flat [23j. The curvature of 
the field-space manifold, in turn, can induce additional interactions among the matter fields, 
beyond those included in the Jordan-frame potential. Moreover, the scalar fields necessarily 
acquire noncanonical kinetic terms in the Einstein frame. These new features can have a 
dramatic impact on the behavior of the fields during inflation, and hence on the primordial 
power spectrum. 

Chief among the multifield effects for producing new features in the primordial power 
spectrum is the ability of fields' trajectories to turn in field-space as the system evolves. 
Such turns are not possible in single-field models, which include only a single direction of 
field-space. In the case of multiple fields, special features in the effective potential, such as 
ridges or bumps, can focus the background fields' trajectories through field space or make 
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them diverge. When neighboring trajectories diverge, primordial bispectra can be amphfied 
to sufficient magnitude that they should be detectable in the CMB [T2HI51 [21H3n]- 

To date, features like ridges in the effective potential have been studied for the most part 
phenomenologically rather than being strongly motivated by fundamental physics. Here we 
demonstrate that ridges arise naturally in the Einstein-frame effective potential for mod- 
els that incorporate multiple fields with nonminimal couplings. Likewise, as noted above, 
models with multiple nonminimally coupled scalar fields necessarily include noncanonical 
kinetic terms in the Einstein frame, stemming from the curvature of the field-space mani- 
fold. Both the bumpy features in the potential and the nonzero curvature of the field-space 
manifold routinely cause the fields' evolution to depart from slow-roll for some duration of 
their evolution during inflation. 

Recent analyses of primordial non-Gaussianities have emphasized two distinct types of 
fine-tuning needed to produce observable bispectra: fine-tuning the shape of the effective po- 
tential to include features like ridges; and separately fine-tuning the fields' initial conditions 
so that the fields begin at or near the top of these ridges [27H30] . Here we show that the ffist 
of these types of fine-tuning is obviated for multifield models with nonminimal couplings; 
such features of the potential are generic. The second type of fine-tuning, however, is still 
required: even in the presence of ridges and bumps, the fields' initial conditions must be 
fine-tuned in order to produce measurable non-Gaussianities. 

In Section II we examine the evolution of the fields in the Einstein frame and emphasize 
the ubiquity of features such as ridges that could make the fields' trajectories diverge in 
field space. Section HI introduces our covariant, multifield formalism for studying the evo- 
lution of background fields and linearized perturbations on the curved field-space manifold. 
In Section IV we analyze adiabiatic and entropy perturbations and quantify their coupling 
using a covariant version of the familiar transfer-function formalism flU\ \T2\ [31] . In Section 
V we build on recent work [32] - [M] to calculate the primordial bispectrum for multifield 
models, applying it here to models with nonminimal couplings. We find that although the 
nonminimal couplings induce new interactions among the entropy perturbations compared 
to models in which all fields have minimal coupling, the dominant contribution to the bis- 
pectrum remains the familiar local form of Jnl, made suitably covariant to apply to the 
curved field-space manifold. Concluding remarks follow in Section VI. We collect quantities 
relating to the curvature of the field-space manifold in the Appendix. 
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II. EVOLUTION IN THE EINSTEIN FRAME 

We consider TV scalar fields in (3 + 1) spacetime dimensions, with spacetime metric 
signature (— , +, +, +). We work in terms of the reduced Planck mass, Mpi = {SnG)^^/'^ = 
2.43 X 10^® GeV. Greek letters label spacetime indices, /i, z/ = 0,1,2,3; lower-case Latin 
letters label spatial indices, i,j = 1,2,3; and upper-case Latin letters label field-space 
indices, /, J = 1, 2, ...,Af. 

In the Jordan frame, the scalar fields' nonminimal couplings to the spacetime Ricci cur- 
vature scalar remain explicit in the action. We denote quantities in the Jordan frame with 
a tilde, such as the spacetime metric, g^u{x). The action for Af scalar fields in the Jordan 
frame may be written 



•5* Jordan I d X 



(1) 



where is the nonminimal coupling function and V"(0^) is the potential for the scalar 

fields in the Jordan frame. We have included the possibility that the scalar fields in the Jor- 
dan frame have noncanonical kinetic terms, parameterized by coefficients Gij{(f)^). Canon- 
ical kinetic terms correspond to Qu = 6jj. 

We next perform a conformal transformation to work in the Einstein frame, in which 
the gravitational portion of the action assumes Einstein-Hilbert form. We define a rescaled 
spacetime metric tensor, g^u{x), via the relation, 

g^^{x) = Vl?{x)g^y{x), (2) 

where the conformal factor is related to the nonminimal coupling function as 
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^'(^) = T7y/(</''(^))- (3) 



Ml, 



2 2 

The potential in the Einstein frame is scaled by the conformal factor 

1 Mt 



Eq. ([T]) then takes the form [23| 

'S'sinstein = d X\/ —g 



(4) 
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Ml 



(6) 



The coefficients of the noncanonical kinetic terms in the Einstein frame depend on the 
nonminimal couphng function, f{(f)'), and its derivatives, and are given by 

3 

- Cr 1 -I 

2/(0^ 
where fj = df /dcj)^. 

As demonstrated in ^3], the nonminimal couphngs induce a field-space manifold in the 
Einstein frame, associated with the metric Qij{(j)^) in Eq. ([6|, which is not conformal to flat 
for models in which multiple scalar fields have nonminimal couplings in the Jordan frame. 
Thus there does not exist any combination of conformal transformation plus field-rescalings 
that can bring the induced metric into the form Qu = Sjj. In other words, multifield models 
with nonminimal couplings necessarily include noncanonical kinetic terms in the Einstein 
frame, even if the fields have canonical kinetic terms in the Jordan frame, Qu = 6ij. When 
analyzing multifield inflation with nonminimal couplings, we therefore must work either with 
a noncanonical gravitational sector or with noncanonical kinetic terms. Here we adopt the 
latter. Because there is no way to avoid noncanonical kinetic terms in the Einstein frame 
in such models, we do not rescale the fields. For the remainder of the paper, we restrict 
attention to models with canonical kinetic terms in the Jordan frame, Qu = 6u, in which 
the curvature of the field-space manifold in the Einstein frame depends solely upon f{(f)^) 
and its derivatives. 

Varying the action of Eq. Q with respect to g^u{x) yields the Einstein field equations, 

Rfiu — -^dfiuR = 'j^T^^, (7) 



where 

~1 



^Qij9"^d^<j)'dp<jy' + V{<jy') 



(8) 



Tf,u = Qijd^(j)^d^(j)-^ - g,,^ 
Varying Eq. Q with respect to (j)^ yields the equation of motion, 

□0^ + gf^-'T'jj.d^cP'd,^'' - Q'^'Vk = 0, (9) 

where = g^'^(j)^.^.^, and r^j^(0^) is the Christoffel symbol for the field-space manifold, 
calculated in terms of Qjj. 

We expand each scalar field to ffist order around its classical background value, 

<p\xn=^\t) + 6<p\xn, (10) 
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and also expand the scalar degrees of freedom of the spacetime metric to first order, per- 
turbing around a spatially flat Friedmann- Robertson- Walker (FRW) metric [TUl [TTl [^2] , 

(„) 

= - (1 + 2A) dt^ + 2a (diB) dx'dt + [(1 - 2ip)5ij + 2didjE] dx'dx^, 

where a{t) is the scale factor. To background order, the 00 and ij components of Eq. ([T]) 
may be combined to yield the usual dynamical equations, 

'1 



H 



2 



3M2 



1 



(12) 



2My ■'^ ^ ' 



where H = d/a is the Hubble parameter, and the field-space metric is evaluated at back- 
ground order, Qu = Qui.^^)- 

Both the curvature of the field-space manifold and the form of the effective potential in 
the Einstein frame depend upon the nonminimal coupling function, /(0^). The requirement 
of renormalizability for scalar matter fields in a (classical) curved background spacetime 
dictates the form of /(0^) [I7H20]: 



' 2 



(13) 



I 

where Mq is some mass-scale that could be distinct from Mpi, and the nonminimal couplings 
^7 are dimensionless constants that need not be equal to each other. If any of the fields 
develop nonzero vacuum expectation values, {(j)^) = , then one may expect Mpj = Mq + 
C,i{v^Y . Here we will assume either that = for each field or that ^f^iv^ <^ Mpi, so 
that Mo ~ Mpi. 

The nonminimal couplings S,i could in principle take any "bare" value. (Conformal cou- 
plings correspond to ^/ = —1/6; we only consider positive couphngs here, ^/ > 0.) Under 
renormalization-group flow the constants vary logarithmically with energy scale. The exact 
form of the (3 functions depends upon details of the matter sector, but for models whose 
content is akin to the Standard Model the (3 functions are positive and the flow of ^/ has 
no fixed point, rising with energy scale without bound [20]. Studies of the flow of ^ in the 
case of Higgs inflation [36J indicate growth of ^ by (9(10^ — 10^) between the electroweak 
symmetry-break scale, A ~ 10^ GeV, and typical inflationary scales, A ~ 10^^ GeV [37] . 
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Hence we anticipate that realistic models will include nonminimal couplings ^ 1 during 
inflation. 

Renormalizable potentials in (3 + 1) spacetime dimensions can include terms up to quartic 
powers of the fields. A potential in the Jordan frame that assumes a generic renormalizable, 
polynomial form such as 

y = ^ E i<t^') m ' (0^) m " (14) 

/ KJ I 

will yield an effective potential in the Einstein frame that is stretched by the conformal 
factor in accord with Eq. (|5|. As the Jth component of (p^ becomes arbitrarily large the 
potential in that direction will become asymptotically flat, 



(no sum on J), unlike the quartic behavior of the potential in the large- field limit in the 
Jordan frame. (The flatness of the effective potential for large field values was one inspiration 
for Higgs inflation |36].) Inflation in such models occurs in a regime of fleld values such that 
ij{^''f > M^j for at least one component, J. As emphasized in j3Hj, for large nonminimal 
couplings, ^ 1, all of inflation therefore may occur for fleld values that satisfy \^'^\ < Mpi, 
unlike the situation for ordinary chaotic inflation with polynomial potentials and minimal 
couplings. 

Although the effective potential in the Einstein frame will asymptote to a constant value 
in any given direction of fleld space, the constants will not, in general, be equal to each 
other. Thus at flnite values of the flelds, the potential will generically develop features, 
such as ridges or bumps, that are absent from the Jordan-frame potential. Because the 
asymptotic values of V{(l)^) in any particular direction are proportional to 1/0, the steepness 
of the ridges depends sharply on the ratios of the nonminimal coupling constants. If some 
explicit symmetry, such as the SU{2) electroweak gauge symmetry obeyed by the Higgs 
multiplet in Higgs inflation (SB], forces all the couplings to be equal — C,i = ^, Tn'j = m^, 
and A/ = gu = A for all I, J — then the ridges in the Einstein-frame potential disappear 
and the potential asymptotes to the same constant value in each direction of fleld space. 
We study the dynamics of such special cases in |38]. For the remainder of this paper, we 
consider models in which the constants are of similar magnitude but not exactly equal to 
each other. 
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FIG. 1: The Einstein-frame effective potential, Eq. (18), for a two-field model. The potential 



shown here corresponds to the coupHngs ixl^4> ~ 0-^' ^xl^4> ~ 0-^' 9/^4> — O-l' 



10-2 A^Mp2 . 



For definiteness, consider a two-field model with a potential in the Jordan frame of the 
form 



X 



2 ?>r . 2 

and nonminimal coupling function given by 

/(0,X) = J[M2 + 



4 



(16) 



(17) 



In the Einstein frame the potential becomes 
See Fig. [1} 

In addition to the ridges shown in Fig. [T| other features of the Einstein-frame potential 
can arise depending on the Jordan-frame couplings. For example, the tops of the ridges 
can develop small indentations, such that the top of a ridge along x ~ becomes a local 
minimum rather than a local maximum. In that case, field trajectories that begin near the 
top of a ridge tend to focus rather than diverge, keeping the amplitude of non-Gaussianities 



very small. For the two-field potential of Eq. (18), we find 



'lx=o 



m: 



(19) 
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For realistic values of the masses that satisfy m|, <^ Mpj, and at early times when 
'-^ ^ Mpi, the top of the ridge along the x ~ direction will remain a local maximum if 



< Kix- (20) 



When the couplings satisfy Eq. (20), the shape of the potential in the vicinity of its ridges 
is similar to that of the product potential, V = m?e~'^'^^ which has been studied in detail 
in [2H1 EO]- Trajectories of the fields that begin near each other close to the top of a ridge 
will diverge as the system evolves; that divergence in trajectories can produce a sizeable 
amplitude for the bispectrum, as we will see below. 

Even potentials with modest ratios of the nonminimal couplings can produce trajectories 
that diverge sharply, as shown in Fig. |2} As we will see in Section V, trajectory 2 of Fig. [2] 
(solid red line) yields a sizeable amplitude for the bispectrum that is consistent with present 
bounds, whereas trajectories 1 and 3 produce negligible non-Gaussianities. We will return 
to the three trajectories of Fig. 2 throughout the paper, as illustrations of the types of field 
dynamics that yield interesting possibilities for the power spectrum. 



Unlike the product potential studied in |281 130] , the potential of Eq. ( 18 ) contains valleys 
in which the system will still inflate. For trajectories 1 (orange dotted line) and 2 (red solid 
line) in Fig. [2| for example, the system begins near x ~ ^-^id rolls off the ridge; because 
\xl^x 7^ 0' valleys in the x direction are also false vacua and hence the system continues 
to inflate as the fields relax toward the global minimum at = x = 0. Near the end of 
inflation, when ^^0^ + ■CxX^ < ^^e fields oscillate around the global minimum of the 
potential, which can drive a period of preheating. See Fig. [3j 

Evolution of the fields like that shown in Fig. |3]is generic for this class of models when the 
fields begin near the top of a ridge, and can produce interesting phenomenological features 
in addition to observable bispectra. For example, the oscillations of around = when 
the system first rolls off the ridge could produce an observable time-dependence of the 
scale factor during inflation, as analyzed in [ID]. The added period of inflation from the 
false vacuum of the x valley could lead to scale-dependent features in the power spectrum 
associated with double-inflation [41J. 

In the class of models we consider here, neighboring trajectories may also diverge if we 
include small but nonzero bare masses for the fields. For example, in Fig. |4] we show the 
evolution of the fields for the same initial conditions as trajectory 3 of Fig. |2] — the black, 
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FIG. 2: Parametric plot of the fields' evolution superimposed on the Einstein-frame potential. 
Trajectories for the fields (j) and x that begin near the top of a ridge will diverge. In this case, the 

couplings of the potential are = 10, ^y. = 10.02, A^^/A^ = 0.5, g/X^p = 1, and m<^ = rriy- = 0. 

(We use a dimensionless time variable, r = Mpi t, so that the Jordan-frame couplings are 
measured in units of A^.) The trajectories shown here each have the initial condition 0(ro) = 3.1 
(in units of Mpi) and different values of x('^o): xi^o) = 1-1 x 10~^ ("trajectory 1," yellow dotted 
line); xi^o) = 1-1 x 10~^ ("trajectory 2," red solid line); and x{to) = 1-1 x 10""^ ("trajectory 3," 



black dashed line). 




T 



FIG. 3: The evolution of the Hubble parameter (black dashed line) and the background fields, 
(P{t) (red solid line) and xi'T) (blue dotted line), for trajectory 2 of Fig. 2. (We use the same 
units as in Fig. 2, and have plotted lOOH so its scale is commensurate with the magnitude of the 
fields.) For these couplings and initial conditions the fields fall off the ridge in the potential at 
r = 2373 or = 66.6 efolds, after which the system inflates for another 4.9 efolds until 

Tend = 2676, yielding iVtotai = 71.5 efolds. 
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FIG. 4: Models with nonzero masses include additional features in the Einstein- frame potential 
which can also cause neighboring field trajectories to diverge. In this case, we superimpose the 
evolution of the fields (f) and x on the Einstein-frame potential. The parameters shown here are 
identical to those in Fig. [2] but with = 0.075 M^^ and = 0.0025 M^j rather than 0. 
The initial conditions match those of trajectory 3 of Fig. [2| (p{To) = 3.1 and xi^o) = 1-1 ^ in 

units of Mpi. 



dashed curve that barely deviates from the middle of the ridge. The evolution shown in Fig. 
[2] was for the case = = 0. If, instead, we include nonzero masses, then the curvature 
of the effective potential at small field values becomes different from the zero-mass case. In 
particular, for positive, real values of the masses, the ridges develop features that push the 
fields off to one side, recreating behavior akin to what we found in trajectories 1 and 2 of 
Fig. |2} 

Because the field-space manifold is curved, the fields' trajectories will turn even in the 
absence of tree-level couplings from the Jordan-frame potential: the fields' geodesic motion 
alone is nontrivial. The Ricci scalar for the field-space manifold in the two-field case is given 



in Eq. (117). In Fig. [5jwe plot the fields' motion in the curved manifold for the case when 
x) = ^(0) x) = 0. The curvature of the manifold is negligible at large field values but 
grows sharply near ~ x ~ 0. 



Given the nonvanishing curvature of the field-space manifold, we must study the evolution 
of the fields and their perturbations with a covariant formalism, to which we now turn. 
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FIG. 5: Parametric plot of the evolution of the fields and x superimposed on the Ricci 
curvature scalar for the field-space manifold, 7^, in the absence of a Jordan-frame potential. The 
fields' geodesic motion is nontrivial because of the nonvanishing curvature. Shown here is the 
case = 10, = 10.02, (/)(to) = 0.75, x(to) = 0.01, 0'(ro) = -0.01, and x'(to) = 0.005. 



III. COVARIANT FORMALISM 



A gauge-invariant formalism for studying perturbations in multifield models in the Jordan 
frame was developed in ^21 113] . In this paper we work in the Einstein frame, following the 
approach established in [BHSl ESI EHl E2H31]- Our approach is especially indebted to the 
geometric formulation of [30]. In [30], the authors introduce a particular tetrad construciton 
with which to label the field-space manifold locally, which they dub the "kinematical basis." 
The adoption of the kinematical basis simplifies certain expressions and highlights features 
of physical interest in the primordial power spectrum, but it does so at the expense of 
obscuring the relationship between observable quantities and the fields that appear in the 
original Lagrangian, in terms of which any given model is specified. Rather than adopt the 
kinematical basis here, we develop a covariant approach in terms of a single coordinate chart 
that covers the entire field manifold. This offers greater insight into the global structure 
of the manifold, as illustrated in Fig. [5j We also keep coordinate labels explicit, which 
facilitates application of our formalism to the original basis of fields, that appears in the 



13 



governing Lagrangian. Also unlike [30], we work in terms of cosmic time, t, rather than the 
number of efolds during inflation, A^, because we are interested in applying our formalism 
(in later work) to eras such as preheating, for which is a poor dynamical parameter. 
Because of these formal distinctions from [30] , we briefly introduce our general formalism in 
this section. 

We expand each scalar field to first order around its classical background value, as in 



Eq. (10). The background fields, ^p\t), parameterize classical paths through the curved 
field-space manifold, and hence can be used as coordinate descriptions of the trajectories. 
Just like spacetime coordinates in general relativity, x^, the array ip^ is not a vector in the 
field-space manifold [44J. Infinitesimal displacements, dip^ , do behave as proper vectors, and 
hence so do derivatives of (f^ with respect to an affine parameter such as t. 

For any vector in the field space, , we define a covariant derivative with respect to the 
field-space metric as usual by 

VjA' = djA' + T'jkA''. (21) 

Following [71 [251 EO], we also introduce a covariant derivative with respect to cosmic time 
via the relation 

VtA' = ip^VjA' = A' + T'j^A^if'', (22) 



where overdots denote derivatives with respect to t. The construction of Eq. (22) is essen- 
tially a directional derivative along the trajectory. 

For models with nontrivial field-space manifolds, the tangent space to the manifold at 
one time will not coincide with the tangent space at some later time. Hence the authors of 
[32] l33] introduce a covariant means of handling field fluctuations, which we adopt here. As 



specified in Eq. (10), the value of the physical field at a given location in spacetime, 0^(x^), 
consists of the homogenous background value, f^(t), and some gauge-dependent fluctuations, 
6(j)\x^). The fluctuation 6(f)' represents a finite coordinate displacement from the classical 
trajectory, and hence does not transform covariantly. This motivates a construction of a 
vector Q' to represent the field fluctuations in a covariant manner. The two field values, 
(p' and ip' , may be connected by a geodesic in the field-space manifold parameterized by 
some parameter A, such that = 0) = ip' and 4>'{\ = 1) = ip' + 6(j)' . These boundary 
conditions allow us to identify a unique vector, Q' , that connects the two field values, such 
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that V- 



\V |A=0 



. One may then expand 5(j)^ in a power series in j32| 133] . 



(23) 



where the Christoffel symbols are evaluated at background order in the fields, F^j^ = 
jj^{ip^). To first order in fiuctuations 5(f)^ — ?■ Q^, and hence at linear order we may treat 
the two quantities interchangeably. When we consider higher-order combinations of the 
field fiuctuations below, however, such as the contribution of the three-point function of 
field fiuctuations to the bispectrum, we must work in terms of the vector rather than 



We introduce the gauge-invariant Mukhanov-Sasaki variables for the perturbations pUl 

[III Eg, 

(24) 



Because both and (p^ are vectors in the field-space manifold, is also a vector. The 
Mukhanov-Sasaki variables, , should not be confused with the vector of field fiuctuations, 

. The are gauge- invariant with respect to spacetime gauge transformations up to first 
order in the perturbations, and are constructed from a linear combination of field fiuctuations 
and metric perturbations. The quantity does not incorporate metric perturbations; it is 
constructed from the (gauge-dependent) field fiuctuations and background-order quantities 
such as the field-space Christoffel symbols. At lowest order in perturbations, — )• in 
the spatially fiat gauge. 

Using Eq. (24), Eq. ^ separates into background and first-order expressions. 



VtP' + 3if0^ + g'^'VK = 0, 



and 



V^Q' + 3HVtQ' + 



.2" J 



'pi- 



Q' = 0. 



(25) 



(26) 



The mass-squared matrix appearing in Eq. ( 26 ) is given by 



(27) 



where TZ^lmj is the Riemann tensor for the field-space manifold. All expressions in Eqs. 



(25), (26), and (27) involving Qu, T^jj^, TZ^lmj^ ^ evaluated at background order 
in the fields, ip^ . 
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The system simplifies further if we distinguish between the adiabatic and entropic direc- 
tions in field space [H]. The length of the velocity vector for the background fields is given 
by 

(28) 



Introducing the unit vector, 



a 



(29) 



the background equations, Eqs. (|12|) and (|25|), simplify to 

'1 



H 



1 



3M2 



1 



(30) 



2M2 



and 



where we have defined 



a + 3H& + V„ = 0, 



(31) 



(32) 



The background dynamics of Eqs. (30) and (31) take the form of a single-field model with 



canonical kinetic term, with the exception that V{ip^) in Eqs. (30) and (31) depends on all 
independent fields, ip^ . 

The directions in field space orthogonal to are spanned by 



s = y — a a . 

The quantities and s^'^ obey the useful relations 

aici^ = 1, 



(33) 



Sjj=^-1, 



ajs^-^ = for all J. 



(34) 



Therefore we may use and s^"^ as projection operators to decompose any vector in field 
space into components along the direction and perpendicular to as 



(35) 
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FIG. 6: The slow-roll parameters e (blue dashed line) and \r]aa\ (solid red line) versus A^* for 
trajectory 2 of Fig. 2, where A'^* is the number of efolds before the end of inflation. Note that 
\ricrcr\ temporarily grows significantly larger than 1 after the fields fall off the ridge in the potential 

at around A^^, ~ 5. 



In particular, 5*^ = j'p'^ vanishes identically, 5*^ = 0. Thus all of the dynamics of the 
background fields are captured by the behavior of a and . 



Given the simple structure of the background evolution, Eqs. (30) and (31), we introduce 
slow-roll parameters akin to the single-field case. We define 



H 



3(7^ 



((t2 + 2V) 



and 



where we have defined 



M. 



pi" 



V 



(36) 



(37) 



(38) 



M.J = aiM'j = (VkVjV) , 
M.. = ojg'M'j = a^'a' (VkVjV) . 

The term in Ai^j involving TZ^lmj vanishes when contracted with ai or o"^ due to the first 
Bianchi identity (since the relevant term is already contracted with and hence Aiaa 

is independent of TI^imj- -^^^ trajectory 2 of Fig. |2] (solid red line), we see that slow- roll is 
temporarily violated when the fields roll off the ridge of the potential. See Fig. [6j 

A central quantity of interest is the turn-rate [30], which we denote . The turn-rate is 
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FIG. 7: The turn-rate, co = \u}^\, for the three trajectories of Fig. 2] trajectory 1 (orange dotted 
hne); trajectory 2 (red sohd hne); and trajectory 3 (black dashed hne). The rapid oscillations in 

id correspond to oscillations of (j) in the lower false vacuum of the x valley. For trajectory 1, uj 
peaks at = 34.5 efolds before the end of inflation; for trajectory 2, u peaks at iV* = 4.9 efolds 
before the end of inflation; and for trajectory 3, cj remains much smaller than 1 for the duration 

of inflation. 



given by the (covariant) rate of change of the unit vector, , 



1 

a 



(39) 



where the last expression follows upon using the equations of motion, Eqs. (25) and (31). 
Because oc s^^ , we have 

oj^ai = 0. (40) 



Using Eqs. (33) and (39), we also find 



(41) 



For evolution of the fields like that shown in Fig. |2| the turn-rate peaks when the fields roll 



off the ridge; see Fig. 39 



We may decompose the perturbations along directions parallel to and perpendicular to 



5s^ = s^jQ^. 



(42) 
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Note that Ss^ may be defined either in terms of the field fluctuations or the Mukhanov- 
Sasaki variables, since s^jScj)"^ = s^jQ"^- Though 6s^ is a vector in field-space with J\f 
components, only A/" — 1 of these components are linearly independent. We will isolate 
particular components of interest in Section IV. 

Taking a Fourier transform, such that for any function F{t, x*) we have a'^{t)did'^F{t, x*) = 



—k'^Fk{t) where k is the comoving wavenumber, Eq. (26) separates into two equations of 
motion (we suppress the label k on Fourier modes). 



k^ 



d 



Ml^a^ dt 




(43) 



and 



Ct \ (T 



-2uj' 



H a 



(44) 



Although the effective mass of the adiabatic perturbations, m^g = A^o-o- — i^^, is independent 

of 

LMj^ the curvature of the field-space manifold introduces couplings among components 
of the entropy perturbations, 5s^ , by means of the M.^ j term in Eq. (44). The quantities 



and {(jjjds'^) are scalars in field space, so the covariant time derivatives in Eq. (43) reduce 
to ordinary time derivatives. 



From Eqs. (43) and (44), it is clear that the adiabatic and entropy perturbations decouple 



if the turn-rate vanishes, uj^ = 0. Moreover, Eq. (43) for Q^r is identical in form to that of 
a single-field model (with m^g = Ai^a — <-^^), but with a nonzero source term that depends 
on the combination ujj5s^ . Even in the presence of significant entropy perturbations, , 
the power spectrum for adiabatic perturbations will be devoid of features such as non- 
Gaussianities unless the turn-rate is nonzero, ^ 0. 



IV. ADIABATIC AND ENTROPY PERTURBATIONS 



In Section III we identified the vector of entropy perturbations, 5s^ , which includes M —1 
physically independent degrees of freedom. As we will see in this section, these Af — 1 
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physical components may be further clarified by introducing a particular set of unit vectors 
and projection operators in addition to and s^'^ . With them we may identify components 
of 5s^ of particular physical interest. 

We denote the gauge-invariant curvature perturbation as IZc, not to be confused with the 
Ricci scalar for the field-space manifold, IZ. The perturbation Tic is defined as (TUl [HI E2] 

nc = ^- r^M, (45) 

where p and p are the background-order energy density and pressure for the fiuid filling the 
FRW spacetime, and 5q is the energy-density fiux of the perturbed fiuid, T° = di5q. Given 
Eq. g, we find 

P-]^^^V. 

P=\^'- V. (46) 



and hence, upon using Eqs. (24) and (42) 



7^c = ^ + = -Q.. (47) 

a (7 



We thus find that TZc oc Qa, and that the righthand side of Eq. (44) is proportional to Tic- 
Recall that these expressions hold to first order in fiuctuations, for which 50^ — > . 

In the presence of entropy perturbations, the gauge-invariant curvature perturbation need 
not remain conserved. Tic 7^ 0. In particular, the nonadiabatic pressure perturbation is given 
by [mill] 

5pnad = 6p- %5p = -^rrfr^p^ + 2a iu}j5s'^) , (48) 
p 3Ha ^ ' 

where 8pm = Sp—diHSq is the gauge-invariant comoving density perturbation. The perturbed 

Einstein field equations (to linear order) require [TOl |Tl] 

5p^ = -2Mp2 (49) 

where \1/ is the gauge-invariant Bardeen potential (TUl HH 

= i) + a'H (e . (50) 
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Therefore in the long- wavelength limit, for k <^ aH, the only source of nonadiabatic pressure 
comes from the entropy perturbations, 6s^ . Using the usual relations (TUl [H] among the 
gauge-invariant quantities TZc and ( = —ip + {H/p)Sp, we find 



Ha a 



(51) 



Thus even for modes with k ^ aH, TZc will not be conserved in the presence of entropy 
perturbations if the turn-rate is nonzero, uj^ 7^ 0. 



Eqs. (43) and (51) indicate that a particular component of the vector 5s^ is of special 

physical relevance: the combination {ujjds'') serves as the source for Q^j and hence for TZc- 

Akin to the "kinematical basis" of [30j, we may therefore introduce a new unit vector that 

points in the direction of the turn-rate, uj\ together with a new projection operator that 

picks out the subspace perpendicular to both and uj^: 

J 

(52) 



s = 



IJ 



where u = \u)^\ is the magnitude of the turn-rate vector. Using the relations in Eq. (34), 



the definitions in Eq. (52) imply 



s'rs' + 7", 

Af-2, 



(53) 



We then find 



ajs^ = an^'^ = sa" = for all J. 



IJ 



where 



= -M^Ki 



IK 



and hence, from Eq. (53), 



a/U^ = siH^ = 0. 
The vector of entropy perturbations, 6s^ , may then be written as 

5/ = s^Qs + B\ 



(54) 

(55) 
(56) 

(57) 
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where 



Qs = sjQ^, 



(58) 



The quantity that sources and IZc is now easily identified as the scalar, ojjds'^ = uQs, 
which corresponds to just one component of the vector 6s^ . 



Making use of Eqs. (30), (47), and (51), the equation of motion for in Eq. (44) 
separates into 



ru2 



Qs + SHQs + 



: 4Mp2 -^^ - Vt (UjB') - UjVtB' - MsjB' - 3H (UjB') 



(59) 



and 



V^B^ + [3H5[j + 2 {a^uj - s^Uj)] VtB" 



+ 



2U'Qs - i'^MasQs + {SHU' + VtU') Qs 



B' 



(60) 



In analogy to (38), we have introduced the projections 

Msj = siM'j, 



(61) 



Note, however, that unlike A^o-j, the term in M.' j proportional to Tl' lj^.^j does not vanish 
upon contracting with si or . Hence the Riemann-tensor term in M.' j induces interactions 
among the components of 5s' . 

For models with A/" > 3 scalar fields, we may introduce additional unit vectors and 
projection operators with which to characterize components of B' . The next in the series 
are given by 



U 



n 



(62) 



Repeating steps as before, we find 



Vtii' = Us' + t' 



u'r' - t'u\ 



(63) 
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where 



We then have 



in terms of 



n 



MsK + -cr^cTL [VaM'-k) 



JK 



(64) 



(65) 



Qu = ujQ 



(66) 



This decomposition reproduces the structure in the "kinematical basis" [30] but can be 
apphed in any coordinate basis for the field-space manifold: Q„ is sourced by Qg] Qg is 



sourced by Qo- and Qu (though we have used Eq. (51) to substitute the dependence on Q 



for the V^\E' term in Eq. (59)); Qu is sourced by Qs and Q^ = tjQ'^ /\t^\, and so on. 

For our present purposes the two-field model will suffice. The perturbations then consist of 



two scalar degrees of freedom, Qo- and Qs, which obey Eqs. (43) and (59) (with = = 0) 
respectively. The effective mass-squared of the entropy perturbations becomes 



2 - Mss + 3u;^ 



(67) 



If the entropy perturbations are heavy during slow-roll, with fig > 3if/2, then the amplitude 
of long- wavelength modes, with k ^ aH, will fall exponentially: Qs ~ a^^^'^{t) during quasi- 
de Sitter expansion. For trajectories that begin near the top of a ridge, on the other hand, 
the entropy modes will remain light or even tachyonic at early times, since /i^ is related to the 
curvature of the potential in the direction orthogonal to the background fields' trajectory. 
Once the background fields roll off the ridge, the entropy mass immediately grows very large, 
suppressing further growth in the amplitude of Qs- See Fig. [8j 

The perturbations in the adiabatic direction are proportional to the gauge-invariant cur- 



vature perturbation, as derived in Eq. (47). Following the usual convention [31], we may 
define a normalized entropy perturbation as 



S = -Qs. 
a 



(68) 



In the long- wavelength limit, the coupled perturbations obey general relations of the form 
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FIG. 8: The effective mass-squared of the entropy perturbations relative to the Hubble scale, 
(fis/H)'^, for the trajectories shown in Fig. [2] trajectory 1 (orange dotted line); trajectory 2 (red 
solid line); and trajectory 3 (black dashed line). For all three trajectories, //^ < while the fields 
remain near the top of the ridge, since /x^ is related to the curvature of the potential in the 
direction orthogonal to the background fields' evolution. The effective mass grows much larger 
than H as soon as the fields roll off the ridge of the potential. 



Uc = aHS + O 
S = (5HS + O 



(69) 



in terms of which we may write the transfer functions as 

Tns{U,t)= [ dt' a{t')H{t')Tss{U,t') 



Tss{t*,t) = exp 



dt' (5{t')H{t') 



(70) 



The transfer functions relate the gauge- invariant perturbations at one time, t^,, to their 
values at some later time, t. We take to be the time when a fiducial scale of interest first 
crossed outside the Hubble radius during inflation, defined by a?'{t^)H'^{t^) = k^. In the 
two-field case, both TZc and S are scalars in field space, and hence a, P, T-jis, and Tss are 



also scalars. Thus there is no time-ordering ambiguity in the integral for Tss in Eq. (70). 



In the two-field case, Eq. (51) becomes 

n = 2ujS + o 



(71) 
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Comparing with Eq. (69), we find 



2u{t) 



(72) 



The variation of the gauge-invariant curvature perturbation is proportional to the turn-rate. 



For S we take the long- wavelength and slow-roll limits of Eq. (59): 



Eq. (69) then yields 



^ 3H^ Ha 



Taking the slow-roll limit of Eq. (31 ) for a, we have 



3H& 



(73) 



(74) 



(75) 



Taking a covariant time derivative of both sides, using the definition of in Eq. (39), and 
introducing the slow-roll parameter 

M,, 



M, 



pr 



V 



(76) 



4c^ 



(77) 



we arrive at 

/3 = -2e - 'nss + v. 

where rj^j^ is defined in Eq. (37). For trajectories that begin near the top of a ridge, rjss will 
be negative at early times (like /if), which can yield (3 > 0. In that case, Tssit*,t) will grow 



If one also has a nonzero turn-rate, u — and hence, from Eq. (72), a nonzero a within the 
integrand for T'jis{t^:,t) — then the growing entropy modes will source the adiabatic mode. 
The power spectrum for the gauge-invariant curvature perturbation is defined by [TT] 

(7^,(kl)7^,(k2)) = (27r)3(5(3)(ki + k2)P7^(A;l), 

where the angular brackets denote a spatial average and Pnik) 
power spectrum is then given by 



(78) 



iT^rP. The dimensionless 



(79) 



and the spectral index is defined as 



ra„ = 1 + 



dink 



(80) 
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Using the transfer functions, we may relate the power spectrum at time t^, to its value at 
some later time, t, as 

Vn{k)=Vn{K)[l + T^sit*,t)], (81) 

where k corresponds to a scale that crossed the Hubble radius at some time t > t^. The 
scale-dependence of the transfer functions becomes [HI [121 EOl E] , 

= -a- fjTus, 



H dU 

c c 



1 8Tss 



H dU 

and hence the spectral index for the power spectrum of the adiabatic fluctuations becomes 

Us = nM + ^ sin (2A) (83) 

where 



cosA = — (84) 



Given Eq. (43) in the limit Uj6s = uQs ^ 1, the spectral index evaluated at t^, matches 



the usual single-field result to lowest order in slow-roll parameters [101 [HI HH] : 

ns{U) = l-Qt{U) + 2r]^„{U). (85) 

Scales of cosmological interest first crossed the Hubble radius between 40 and 60 efolds 
before the end of inflation. In each of the scenarios of Fig. [2] the fields remained near 
the top of the ridge in the potential until fewer than 40 efolds before the end of inflation. 
As indicated in Fig. |9| T-jis remains small between A^^, = 60 and 40 for each of the three 
trajectories, with little sourcing of the adiabatic perturbations by the entropy perturbations. 



This behavior of T-jis is consistent with the behavior of w = aH/2 as shown in Fig. 39 
u (and hence a) remains small until the fields roll off the ridge in the potential. Only in 
the case of trajectory 1, which began least high on the ridge among the trajectories and 
hence fell down the ridge soonest (at N^, = 34.5 efolds before the end of inflation), does T-jis 
become appreciable by A^,,, = 40. In particular, we find T-jis{N4q) = 0.530 for trajectory 1; 
Ttis{N4o) = 0.011 for trajectory 2; and T-jisiN^o) = 0.001 for trajectory 3. 

Fixing the fiducial scale fc* to be that which first crossed the Hubble radius A^* = 60 
efolds before the end of inflation, we find ^^(t*) = 0.967 for each of the three trajectories of 
Fig. [2| in excellent agreement with the observed value = 0.968 ±0.012 [2]. Corrections to 
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FIG. 9: The transfer function T-jis for the three trajectories of Fig. [2j trajectory 1 (orange 
dotted hne); trajectory 2 (red sohd hne); and trajectory 3 (black dashed hne). Trajectories 2 and 
3, which begin nearer the top of the ridge in the potential than trajectory 1, evolve as essentially 
single-field models during early times, before the fields roll off the ridge. 



rig from the scale-dependence of Ttis remain negligible as long as Ttis remains small between 
A^* = 60 and 40. Consequently, we find negligible tilt in the spectral index across the entire 
observational window for trajectories 2 and 3, whereas the spectral index for trajectory 1 
departs appreciably from ns{U) for scales that crossed the Hubble radius near A^* = 40. See 
Fig. [To 



V. PRIMORDIAL BISPECTRUM 

In the usual calculation of primordial bispectra, one often assumes that the field fluctua- 
tions behave as nearly Gaussian around the time t^,, in which case the three-point function 



for the field fluctuations should be negligible. Using the construction of Eq. (23), the 
authors of [321 133] calculated the action up to third order in perturbations and found several 
new contributions to the three-point function for field fluctuations, mediated by the Rie- 
mann tensor for the field space, TZ^jkl- The presence of the new terms is not surprising; 
we have seen that Ti^ jkl induces new interactions among the perturbations even at linear 



order, by means of the mass-squared matrix, A^^j in Eq. (27). Evaluated at time t^,, the 



27 



Us 

0.98 h 



0.96 
0.94 
0.92 
0.9 
0.88 
0.86 



60 



55 



50 



45 



40 



N. 



FIG. 10: The spectral index, n^, versus N^, for the three trajectories of Fig. [2] trajectory 1 
(orange dotted Hne); trajectory 2 (red sohd hne); and trajectory 3 (black dashed line). The 
spectral indices for trajectories 2 and 3 coincide and show no tilt from the value ns{NQQ) = 0.967. 



three-point function for calculated in [33] takes the form 



(86) 



Upon using the definition of in Eq. (29), the background equation of Eq. (30) to relate 
H to 0"^, and the definition of e in Eq. (36), the terms on the righthand side of Eq. (86) 
may be written 



^IJK 



2e 



Mpi 

4Mpiy2i 
3 

2M2 e 



f^(]^^^ k2, ks) + cychc permutations, 
a/^^^'^^^'^ J'bO^i, k2, ks) + cyclic permutations. 



(87) 



'pi 



3 

8M2i e 



o-actbTZ^^^^^^'^'^^ fci^i, k2, kg) + cychc permutations. 



aAO^sT^^^'^^^^'^ fvO^i, k2, kg) + cyclic permutations. 



where Jl^^^'^'^ = Q^^^Vm'TI^^^'^ ■, and /r(kj) are shape-functions in Fourier space that de- 
pend on the particular configuration of triangles formed by the wavevectors kj. Comparable 
to the findings in |26l |27] , each of the contributions to the three-point function for the field 
fluctuations is suppressed by a power of the slow-roll parameter, e <^ 1. 
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The quantity of most interest to us is not the three-point function for the field fluctua- 
tions but the bispectrum for the gauge-invariant curvature perturbation, (, which may be 
parameterized as 



(C(ki)C(k2)C(k3)) = (27r)35(3) (ki + k2 + ks) i?c(ki, k2, ka). 



(88) 



Recall that the two gauge-invariant curvature perturbations, TZc and (, coincide in the long- 
wavelength limit when working to first order in metric perturbations [TOl [11] . In terms of 
, the 6N expansion for ( on super-Hubble scales becomes 



(89) 



where = In |a(tcnd)-f^(i^cnd)/^*| is the number of efolds after a given scale K first crossed 



the Hubble radius until the end of inflation. At t*, Eqs. (86) and (89) yield 



(C(ki)C(k2)C(k3)). = NjNjNK{Q\kr)Q''{k,)Q''{k,)), 

+ ^{VjVjN)NkNl 



(2^ 



;Q^(ki -q)Q^(k2)).(Q-^(q)Q^(k3)). + cyclic perms. 



(90) 



The bottom two lines on the righthand side give rise to the usual form of /atl, made suit- 
ably covariant to reflect Qu ^ Sjj. Adopting the conventional normalization, this term 
contributes PTI - IT^ : 



(C(ki)C(k2)C(k3))/,, = (27r)35(3) (ki + k2 + ks 



hiS hiS hi3 



X 



[kl + kl + kl] 



(91) 



where 



5 N'^N'^VaVbN 



NL 



6 (NjN'^Y 



(92) 



The term on the first line of Eq. (90), proportional to the nonzero three-point function for 



the field fluctuations, yields new contributions to the bispectrum. However, the three-point 
function {Q^ Q"^ Q^)* is contracted with the symmetric object, NjN^jN^k- Hence we must 
consider each term within A^'^^ with care. 
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The term proportional to A^"^^ contributes 



(C(ki)C(k2)C(k3))^ = (27r)-^(5(3)(ki + k, + k. 



41.31,3^3 



X 



pi 



{NjN^jN^k) [a'g''' + a-'g;^^ + a^g''] Uk) 



(27r)'V=^)(ki + k2 + k3 



Ht 



(93) 



4^3l3l3 



pi 



where is shorthand for /_4(ki, k2, ks). Using Eqs. (22), (29), (30), and (36), we have 

1 



a^Ni = -u^^ViN = -VtN = — = =, 



(94) 



and hence 



Ht 



(C(kOC(k2)C(k3))^ = (2vr)=^5(3)(ki + k2 + k3)^7:3T:3p 



Ml 



Uk). (95) 



The term arising from B^"^^ contributes 



Ht 



(C(kOC(k2)C(k3))B = (27r)35(3)(ki + k^ + 

g a^iV,,iV,^iV,^ [7^^■^^^ + 7^^^■^^ + cychc] Uk). 



4MpiV2e 



(96) 



But from the symmetry properties of the Riemann tensor we have 71^-^^^ = Jl^^^J 
_j^AKij ^ and from the first Bianchi identity. 



7^ 



A[KIJ] 



0. 



(97) 



The antisymmetry of the Riemann tensor in its last three indices means that any contraction 
of the form 

OukTZ^^^'' = (98) 
for objects Ojjk that are symmetric in the indices J, J, K. In our case, we have Ojjk = 



NjNjN^K and thus every term in the square brackets of Eq. (96) including the cyclic 



permutations may be put in the form of Eq. (98). We therefore find 



(C(ki)C(k2)C(k3))i5 = 



(99) 



identically. 
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The term arising from C^^^ contributes 



(C(ki)C(k2)C(k3))c = (27r)35(3)(k^ + k2 + k3 



41,31,31,3 



X 



2M2 e 



(100) 



We may identify the nonzero terms within this expression using the Bianchi identities. The 



first Bianchi identity is given in Eq. (97), and the second Bianchi identity may be written 



7^ 



AB[CD;E\ 



0. 



(101) 



Using the (anti) symmetry properties of the Riemann tensor and Eqs. (97) and (101), to- 



gether with the fact that the combinations Ouk = NjN^jN^x and ^Iab = (^a^b are sym- 



metric in their indices, we find the only nonzero term within Eq. ( 100 ) to be 

(C(ki)C(k2)C(k3))c = (27r)35(3)(ki + k2 + ks^ 

2M2 e 



zip 1,3 1,3 



aA(rBNjNjN,Kn'^'''-'''fc{k). 
The final term to consider arises from I)^'^^ . In particular, we have 

(C(ki)C(k2)C(k3))l, = i27T)'5^'\k^ + k2 + k3^ 



ZlJ,3 1,3^3 



AMI ^ 



(102) 



cjAaBNjN^jN^K [7^"^^^^ + Tl'^'^-^^ + cychc] fv{k). 

(103) 



Again we may use TZjjka = T^kaij = —T^akij and Eq. (97) to put the first term in square 



brackets in Eq. ( 103 ) in the form 



OukTZ 



AKIJ:B 







(104) 



for OjjK symmetric. The same occurs for the second term in square brackets in Eq. (103) 
and for all cyclic permutations of /, J, K. Hence we find 



(C(ki)C(k2)C(k3))7. = 



(105) 



identically. 
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Writing out the shape functions, /x(kj), for each term p3], the bispectrum becomes 



(C(ki)C(k2)C(k3))* 

= (27r)3(5(3)(ki + k2 + k3 



2^3 2fi3 7^3 
"'l"'2"'3 



6 



X < --fNL{N^AN^^f [kl + kl + ^3^] 



Ml 



-2MM + ^k2-k3 + cychc 



2Ml 6 



-aAcrBNjNiNKTZ 



IABJ:K 



4 

-/., 
9 * 



k^ — ktK^ + cychc 



(106) 



where 



kt = ki + k2 + ks, 
^ ~ kikj 

i<j 



K 



(107) 



and je — 0.577 is the Euler-Mascheroni constant. Only two new contributions appear 
beyond the usual /tvl term, and the term involving 71^ j^l suppressed by the slow-roll 
parameter e ^ 1. 

For the squeezed, local configuration, in which /ci ~ /c2 — K and /ca ~ 0, the shape 
functions for each term contribute as 



[kl + kl + kfi ~ 2k 



-2MM + ^k2 ■ k3 

kt 2 ' ' 



kt 



, + -kf 
k^ ' 3 J 9 * 



7e 



4, 



ktK. 



cyclic 
cyclic 



(2A^ - 3.65)A;^ 



(108) 



Recall that = In \a{t)H{t)/k^\. If we evaluate the bispectrum around the time t^,, then 
~ 0. If we consider shorter length scales, which cross the Hubble radius at some time 
t > t^:, then ~ (9(10), for example, if we consider scales the cross at the end of the 
window of obsevational interest, around = 40 efolds before the end of inflation. Thus the 



shape factor associated with the 7^^^-Sj;_ft: ^gj,j^ j^q^ (106) could dominate at late times; 



though the factor of A^ in /c(ki) is more than compensated by the coefficient e in front of 
that term. For times in the vicinity of on the other hand, each shape factor contributes 
with a comparable magnitude for the squeezed, local state. 
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We may now compare the time-dependent coefficients of each term in Eq. (106). The 
second term, stemming from A^"^^ , is proportional to (A^^^A^'"^), whereas the f^L term is 
multiphed to the square of that term. For models of interest here, in which the potential 
includes ridges, the gradient term is significant. For each of the three trajectories of Fig. 2, 
for example, (N^aN^^) = 0(10^) across the full range A^* = 60 to A^,,, = 40. The gradient 
increases as the ratio of Cx/^^ increases, and hence the Jnl term will dominate the term 
coming from A^"^^ whenever I/nlI > 10^^. 



For the final term in Eq. (106), involving TZ^^^'^'^, we take advantage of the fact that 



for two-field models the Riemann tensor for the field space may be written 

TZabcd = JC{(p') [QacQbd - QadQbc] , (109) 

where JC{4>^) is the Gaussian curvature. In two dimensions, /C(0^) = \lZ{(t)^), where TZ is 
the Ricci scalar. Since VkQab = Qab-k = and K,{4>^) is a scalar in the field space, the 
covariant derivative of the Riemann tensor is simply proportional to the ordinary (partial) 
derivative of the Gaussian curvature, /C. In particular, we find 

aACTBNjNjNKn'^''''-^'' = - (f'NjNj) (Ar,^/C'^) , (110) 

where s^'^ = Q^"^ — a^a^ is the projection operator for directions orthogonal to the adiabatic 



direction. We calculate fC in Eq. (117). At early times, as the system undergoes slow-roll 



inflation, we have ^</>0^ + ^xX^ ^ ^pi- Fo^' the trajectories as in Fig. 2, moreover, the system 
evolves along a ridge such that ^^0^ ^ ^xX^- that case, we find 

^ ^ imiM\ + ^^^^ + + ^^^^^^^ " ~ ^^^^^ 

and hence /C'^ ~ 0. Thus, in addition to being suppressed by the slow-roll factor, e, the 
contribution to the primordial bispectrum from the Jl^^^J'^ term is negligible in typical 
scenarios of interest, because of the weak variation of the Gaussian curvature of the field- 
space manifold around the times A^,,, = 60 to A^,,, = 40 efolds before the end of inflation. 
This matches the behavior shown in Fig. [5] the field-space manifold is nearly fiat until one 
reaches the vicinity of 0, x ~ 0, near the end of inflation. 

For models of the class we have been studying here, we therefore expect the (covariant 
version of the) usual Jnl term to dominate the primordial bispectrum, at least for the 
squeezed, local configuration. Given the weak dependence of the Gaussian curvature, /C(0^), 
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FIG. 11: The non-Gaussianity parameter, IJnlI, for the three trajectories of Fig. [2} trajectory 1 
(orange dotted hne); trajectory 2 (sohd red hne); and trajectory 3 (black dashed hne). Changing 
the fields' initial conditions by a small amount leads to dramatic changes in the magnitude of the 

primordial bispectrum. 



with (j)^ between A^* = 60 and N^, = 40, we do not expect any significant contributions to 
the running of J'nl with scale to come from the curvature of the field-space manifold, given 
the analysis in [31]. 

For the three trajectories of Fig. [2| we find the middle case, trajectory 2, yields a value 
of /nl of particular interest: I/tvlI = 43.3 for fiducial scales k^: that first crossed the Hubble 
radius A^* = 60 efolds before the end of inflation. Note the strong sensitivity of /nl to the 
fields' initial conditions: varying the initial value of xi'^'o) by just |Ax(ro)| = 10~^ changes 
the fields' evolution substantially — either causing the fields to roll off the hill too early 
(trajectory 1) or not to turn substantially in field space at all (trajectory 3) — both of 



which lead to neghgible values for f^i- See Fig. 11 



VI. CONCLUSIONS 



We have demonstrated that multifield models with nonminimal couplings generically 
produce the conditions required to generate primordial bispectra of observable magnitudes. 
Such models satisfy at least three of the four criteria identified in previous reviews of pri- 
mordial non-Gaussianities [16] , namely, the presence of multiple fields with noncanonical 
kinetic terms whose dynamics temporarily violate slow-roll evolution. 
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Two distinct features are relevant in this class of models: the conformal stretching of the 
effective potential in the Einstein frame, which introduces nontrivial curvature distinct from 
features in the Jordan-frame potential; and nontrivial curvature of the induced manifold for 
the field space in the Einstein frame. So long as the nonminimal couplings are not precisely 
equal to each other, the Einstein-frame potential will include bumps or ridges that will tend 
to cause neighboring trajectories of the fields to diverge over the course of inflation. Such 
features of the potential are generic to this class of models, and hence are strongly motivated 
by fundamental physics. 

We have found that the curvature of the potential dominates the effects of interest at 
early and intermediate stages of inflation, whereas the curvature of the field-space manifold 
becomes important near the end of inflation (and hence during preheating). The generic na- 
ture of the ridges in the Einstein-frame potential removes one of the kinds of fine-tuning that 
have been emphasized in recent studies of non-Gaussianities in multifield models, namely, 
the need to introduce potentials of particular shapes [13 EZl ESI ED]- (We are presently 
performing an extensive sweep of parameter space to investigate how Jnl behaves as one 
varies the couplings A/, and m/. This will help determine regions of parameter space 
consistent with current observations.) On the other hand, much as in [151 EZl EHl |30], we 
find a strong sensitivity of the magnitude of the bispectrum to the fields' initial conditions. 
Thus the production during inflation of bispectra with magnitude IJnlI ~ O{50) requires 
fine-tuning of initial conditions such that the fields begin at or near the top of a ridge in the 
potential. 

A subtle question that deserves further study is whether the formalism and results derived 
in this paper show any dependence on frame. Although we have developed a formalism 
that is gauge-invariant with respect to spacetime gauge transformations, and covariant with 
respect to the curvature of the field-space manifold, we have applied the formalism only 
within the Einstein frame. The authors of [13] recently demonstrated that gauge-invariant 
quantities such as the curvature perturbation, (, can behave differently in the Jordan and 
Einstein frames for multifield models with nonminimal couplings. The question of possible 
frame-dependence of the analysis presented here remains under study. 
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APPENDIX A: FIELD-SPACE METRIC AND RELATED QUANTITIES 



Given in Eq. (17) for a two-field model, the field-space metric in the Einstein 

frame, Eq. ([6]), takes the form 



2/ 

V 
V 



1 + 



2a2 



f 



1 + 



/ 



The components of the inverse metric are 

g4'<l' - 

g<t>x — gx<f> - 
gxx = 

where we have defined the convenient combination 

The Christoffel symbols for our field space take the form 



2/ \ 




K) 




C 




) 


^Uix^t^x 








2/^ 


"2/ + 6^202 




c 



x<f> 



r 



<Px 



XX 



C 

^xX 
2/ ' 

C 
C 



f 



<t>x 



x4> 



XX 



2/' 

Ui + QQx ixx 



c 



f 



(112) 



(113) 



(114) 



(115) 



For two-dimensional manifolds we may always write the Riemann tensor in the form 

TIabcd = K^{4>^) [QacQbd - GadGbc] , (116) 



36 



where /C(0^) is the Gaussian curvature. In two dimensions, /C(0^) = ^TZ{(j)^), where 'R-{(p^) 



is the Ricci scalar. Given the field-space metric of Eq. (112), we find 



7^(0^) = 2/C(0O = [(1 + 6e^)(l + 6Q{4f) - C'] . (117) 
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